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Abstract. Solving a system of polynomial equations is a ubiquitous problem in the 
applications of mathematics. Until recently, it has been hopeless to find explicit solutions 
to such systems, and mathematics has instead developed deep and powerful theories 
about the solutions to polynomial equations. Enumerative Geometry is concerned with 
counting the number of solutions when the polynomials come from a geometric situation 
and Intersection Theory gives methods to accomplish the enumeration. 

We use Macaulay 2 to investigate some problems from enumerative geometry, illus- 
trating some applications of symbolic computation to this important problem of solving 
systems of polynomial equations. Besides enumerating solutions to the resulting polyno- 
mial systems, which include overdetermined, deficient, and improper systems, we address 
the important question of real solutions to these geometric problems. 

The text contains evaluated Macaulay 2 code to illuminate the discussion. This is in- 
tended as a chapter in a book on applications of Macaulay 2 to problems in mathematics. 
While this chapter is largely expository, the results in the last section concerning lines 
tangent to quadrics are new. 



1. Introduction 

A basic question to ask about a system of polynomial equations is its number of solu- 
tions. For this, the fundamental result is the following Bezout Theorem. 

Theorem 1.1.1. The number of isolated solutions to a system of polynomial equations 

/l(Xi, ... ,Xn) = f2{Xi, ... ,Xn) = ■■■ = fn{Xl, ■■■ , X„) = 

is hounded by did2 ■ ■ ■ dn, where di := deg/j. If the polynomials are generic, then this 
bound is attained for solutions in an algebraically closed field. 

Here, isolated is taken with respect to the algebraic closure. This Bezout Theorem is a 



consequence of the refined Bezout Theorem of Fulton and MacPherson |jTT], §1.23]. 

A system of polynomial equations with fewer than this degree bound or Bezout number 
of solutions is called deficient, and there are well-defined classes of deficient systems that 
satisfy other bounds. For example, fewer monomials lead to fewer solutions, for which 
polyhedral bounds on the number of solutions are often tighter (and no weaker than) 
the Bezout number, which applies when all monomials are present. When the polynomi- 
als come from geometry, determining the number of solutions is the central problem in 
enumerative geometry. 
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Symbolic computation can help compute the solutions to a system of equations that has 
only isolated solutions. In this case, the polynomials generate a zero-dimensional ideal 
/. The degree of / (dim^ A;[X]/J), which is the number of standard monomials in any 
term order gives an upper bound on the number of solutions, which is attained when I is 
radical. 

Example 1.1.2. We illustrate this discussion with an example. Let /i, /2, /s, and fi be 
random quadratic polynomials in the ring Fioi[?/ii, yi2, 2/21, 1/22] ■ 

11 : R = ZZ/101[yll, yl2, y21, y22] ; 

12 : PolynomialSystem = apply(1..4, i -> 

random(0, R) + randomd, R) + rcindoin(2, R) ) ; 

The ideal they generate has dimension and degree 16 = 2^, which is the Bezout number. 

13 : I = ideal PolynomialSystem; 

03 : Ideal of R 

14 : dim I, degree 1 

04 = (0, 16) 

04 : Sequence 

If we restrict the monomials which appear in the fi to be among 

1, 1/11, Z/12, 1/21, Z/22, l/iiZ/22, and 1/12Z/21, 
then the ideal they generate again has dimension 0, but its degree is now 4. 

15 : J = ideal (random(R~4, R~7) * transpose ( 

matrix{{l, yll, yl2, y21, y22, yll*y22, yl2*y21}})); 

05 : Ideal of R 

16 : dim J, degree J 

06 = (0, 4) 

06 : Sequence 

If we further require that the coefficients of the quadratic terms sum to zero, then the 
ideal they generate now has degree 2. 

17 : K = ideal (random(R~4, R~6) * transpose ( 

matrix{{l, yll, yl2, y21, y22, yll*y22 - yl2*y21»)); 

07 : Ideal of R 

18 : dim K, degree K 

08 = (0, 2) 
08 : Sequence 

In Example |4.4.2| , we shall see how this last specialization is geometrically meaningful. 

For us, enumerative geometry is concerned with enumerating geometric figures of some 
kind having specified positions with respect to general Gxed hgures. That is, counting 
the solutions to a geometrically meaningful system of polynomial equations. We use 
Macaulay 2 to investigate some enumerative geometric problems from this point of view. 
The problem of enumeration will be solved by computing the degree of the (0- dimensional) 
ideal generated by the polynomials. 
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2. Solving systems of polynomials 

We briefly discuss some aspects of solving systems of polynomial equations. For a more 
complete survey, see the relevant chapters in f^, 0]- 

Given an ideal / in a polynomial ring k[X], set V(/) := SpecA;[X]//. When / is 
generated by the polynomials /i, . . . , Jn, gives the set of solutions in affine space to 
the system 

(1) fi{x) = ■■■ = MX) = 

a geometric structure. These solutions are the roots of the ideal I. The degree of a zero- 
dimensional ideal / provides an algebraic count of its roots. The degree of its radical 
counts roots in the algebraic closure, ignoring multiplicities. 

2.1. Excess intersection. Sometimes, only a proper (open) subset of affine space is 
geometrically meaningful, and we want to count only the meaningful roots of I. Often 
the roots V(/) has positive dimensional components that lie in the complement of the 
meaningful subset. One way to treat this situation of excess or improper intersection is 
to saturate J by a polynomial / vanishing on the extraneous roots. This has the effect of 
working in the coordinate ring of the complement of V(/) p|, Exer. 2.3]. 

Example 2.2.1. We illustrate this with an example. Consider the following ideal in 

¥r[x,y]. 

i9 : R = ZZ/7[y, x, MoiiomialDrder=>Lex] ; 

110 : I = ideal (y~3*x~2 + 2*y"2*x + 3*x*y, 3*y~2 + x*y - 3*y) ; 

010 : Ideal of R 

Since the generators have greatest common factor y, I defines finitely many points together 
with the line y = 0. Saturate / by the variable y to obtain the ideal J of isolated roots. 

111 : J = saturated, ideal (y)) 

4 3 2 

011 = ideal (x + x + 3x + 3x, y - 2x - 1) 

011 : Ideal of R 

The first polynomial factors completely in F7[a:], 

112 : factor(J_0) 

012 = (x - 2)(x + l)(x + 2)(x)(l) 
ol2 : Product 

and so the isolated roots of / are (0, 1), (2, 5), (5, 4), and (6, 6). 

Here, the extraneous roots came from a common factor in both equations. A less trivial 



example of this phenomenon will be seen in Section 5.2 



2.2. Elimination, rationality, and solving. Elimination theory can be used to study 
the roots of a zero-dimensional ideal / C k[X]. A polynomial h G k[X] defines a map 
k[y] — k[X] (by y ^ h) and a corresponding projection h: Speck[X] -» A^. The 
generator g{y) G k[y] of the kernel of the map k[y] —>■ k[X]/I is called an eliminant and 
it has the property that V{g) = h{V{I)). When h is a coordinate function Xj, we may 
consider the eliminant to be in the polynomial ring k[xi], and we have {g{xi)) = Jfl k[xi\. 
The most important result concerning eliminants is the Shape Lemma 0. 
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Shape Lemma. Suppose h is a linear polynomial and g is the corresponding eliminant of 
a zero- dimensional ideal I C k[X] with deg(/) = (\.eg{g). Then the roots of I are defined 
in the splitting field of g and I is radical if and only if g is square-free. 

Suppose further that h = Xi so that g = g{xi). Then, in the lexicographic term order 
with Xi < X2 < ■ ■ ■ < Xn, I has a Grobner basis of the form: 

(2) g{xi), X2-g2{xi), Xn-gn{xi), 

where deg{g) > deg{gi) for i = 2, . . . ,n. 

When k is infinite and / is radical, an ehminant g given by a generic hnear polynomial h 
will satisfy deg{g) = deg(J). Enumerative geometry counts solutions when the fixed figures 
are generic. We are similarly concerned with the generic situation of deg{g) = deg(J). 
In this case, eliminants provide a useful computational device to study further questions 
about the roots of /. For instance, the Shape Lemma holds for the ideal of Example |2.2.1 . 
Its eliminant, which is the polynomial J_0, factors completely over the ground field F7, 
so all four solutions are defined in F7. In Section 4.3, we will use eliminants in another 
way, to show that an ideal is radical. 

Given a polynomial h in a zero- dimensional ring k[X]/I, the procedure eliminant (h, 
k[y]) finds a linear relation modulo / among the powers 1, /;,, /i^, . . . ,h'^ of h with d 
minimal and returns this as a polynomial in k[y]. This procedure is included in the 
Macaulay 2 package realroots .m2. 

113 : load "realroots .m2" 

114 : code eliminant 

ol4 = — code for eliminant: 
— realroots .m2 : 65-81 
eliminant = (h, C) -> ( 
Z := C_0; 
A := ring h; 
assert ( dim A == ); 
F := coef f icientRing A; 
assert ( isField F ); 
assertC F == coef f icientRing C ); 
B := basis A; 
d := numgens source B; 
M := f old((M, i) -> M I I 

substitute (contract (B, h~(i+l)), F) , 
substitute (contract (B, 1_A) , F) , 
flatten subsets(d, d)); 
N := ((ker transpose M) ). generators ; 
P := matrix {toList apply(0..d, i -> Z~i)} * N; 
(flatten entries (P)) _0 

) 

ol4 : Net 

Here, M is a matrix whose rows are the normal forms of the powers 1, h, h'^, . . . , h'^ of 
h, for d the degree of the ideal. The columns of the kernel N of transpose M are a basis 
of the linear relations among these powers. The matrix P converts these relations into 
polynomials. Since N is in column echelon form, the initial entry of P is the relation of 
minimal degree. (This method is often faster than naively computing the kernel of the 
map k[Z] A given hj Z h, which is implemented by eliminantNaive (h, Z). 

Suppose we have an eliminant g{xi) of a zero-dimensional ideal / C k[X] with deg{g) = 
deg(/), and we have computed the lexicographic Grobner basis (|]). Then the roots of / 
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are 

(3) {{^i,92{^i),... ,9nm\g{^i) = 0}. 

Suppose now that k = Q and we seek floating point approximations for tlie (complex) 
roots of /. Following this method, we first compute floating point solutions to g{C,) = 0, 
which give all the a;i-coordinates of the roots of /, and then use to find the other 
coordinates. The difficulty here is that enough precision may be lost in evaluating gi{^i) 
so that the result is a poor approximation for the other components ^j. 

2.3. Solving with linear algebra. We describe another method based upon numerical 
linear algebra. When / C k[X] is zero-dimensional, A = k[X]/I is a finite-dimensional 
fc-vector space, and any Grobner basis for / gives an efficient algorithm to compute ring 
operations using linear algebra. In particular, multiplication hj h ^ A is a. linear trans- 
formation rrih : A ^ A and the command regularRep(h) from realroots .m2 gives the 
matrix of rrih in terms of the standard basis of A. 

115 : code regularRep 

ol5 = — code for regularRep: 

— realroots. m2: 97-102 
regularRep = f -> ( 

assert ( dim ring f == ) ; 

b : = basis ring f ; 

k := coef f icientRing ring f ; 

substitute (contract (trcinspose b, f *b) , k) 

015 : Net 

Since the action of A on itself is faithful, the minimal polynomial of is the eliminant 
corresponding to h. The procedure charPoly(h, Z) in realroots. m2 computes the 
characteristic polynomial det{Z ■ Id — rrih) of h. 

116 : code charPoly 

016 = — code for charPoly: 

— realroots. m2: 108-116 
CharPoly = (h, Z) -> ( 

A := ring h; 

F := coef f icientRing A; 

S := F[Z] ; 

Z = value Z ; 

mh := regularRep (h) ** S; 

Idz := S_0 * id_ (S~ (numgens source mh)); 

detCldz - mh) 

) 

ol6 : Net 

When this is the minimal polynomial (the situation of the Shape Lemma), this procedure 
often computes the eliminant faster than does eliminant, and for systems of moderate 
degree, much faster than naively computing the kernel of the map k[Z] ^ A given by 
Z \—>- h. 

The eigenvalues and eigenvectors of rrih give another algorithm for finding the roots of 
/. The engine for this is the following result. 

Stickelberger's Theorem. Let h E A and he as above. Then there is a one-to-one 
correspondence between eigenvectors of rrih and roots ^ of I, the eigenvalue of rrih on 
is the value /i(^) of h at ^, and the multiplicity of this eigenvalue (on the eigenvector 
is the multiplicity of the root ^. 
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Since the linear transformations ruh for h & A commute, the eigenvectors are common 
to all irih. Thus we may compute the roots of a zero-dimensional ideal I C k[X] by first 
computing floating-point approximations to the eigenvectors of rrir^^. Then the root 
^ = (^1, . . . ,^„) of / corresponding to the eigenvector has ith coordinate satisfying 

(4) m^^ ■ = ii-w^. 

An advantage of this method is that we may use structured numerical linear algebra after 
the matrices m^;- are precomputed using exact arithmetic. (These matrices are typically 
sparse and have additional structures which may be exploited.) Also, the coordinates 
are linear functions of the floating point entries of v^, which affords greater precision than 
the non-linear evaluations gi{C,i) in the method based upon elimination. While in principle 
only one of the deg(J) components of the vectors in need be computed, averaging the 
results from all components can improve precision. 



2.4. Real Roots. Determining the real roots of a polynomial system is a challenging 
problem with real world applications. When the polynomials come from geometry, this 
is the main problem of real enumerative geometry. Suppose k (Z M. and / C k[X] is 
zero- dimensional. If g is an eliminant of k[X]/I with deg{g) = deg(J), then the real roots 
of g are in 1-1 correspondence with the real roots of I. Since there are effective methods 
for counting the real roots of a univariate polynomial, eliminants give a naive, but useful 
method for determining the number of real roots to a polynomial system. (For some 



applications of this technique in mathematics, see |]T9, 

The classical symbolic method of Sturm, based upon Sturm sequences, counts the 
number of real roots of a univariate polynomial in an interval. When applied to an 
eliminant satisfying the Shape Lemma, this method counts the number of real roots of 
the ideal. This is implemented in Macaulay 2 via the command SturmSequence(f ) of 
realroots .m2 

il7 : code SturmSequence 

ol7 = — code for SturmSequence: 
— realroots. m2: 120-134 
SturmSequence = f -> ( 

assert ( isPolynomialRing ring f ); 

assert ( numgens ring f === 1 ) ; 

R : = ring f ; 

assert ( char R == ) ; 

X := R_0; 

n := first degree f; 

c := new MutableList from toList (0 .. n) ; 
if n >= then ( 
c#0 = f; 

if n >= 1 then ( 

c#l = diff (x,f ) ; 

scan(2 . . n, i -> c#i = - c#(i-2) 7, c#(i-l)); 

)); 

toList c) 

ol7 : Net 

The last few lines of SturmSequence construct the Sturm sequence of the univariate 
argument /: This is {fo, fi, f2, ■ ■ ■) where /o = /, /i = /', and for i > 1, /j is the normal 
form reduction of —fi-2 modulo fi-i- Given any real number x, the variation of / at 
X is the number of changes in sign of the sequence {fo{x), /i(x), f2{x), . . .) obtained by 
evaluating the Sturm sequence of / at x. Then the number of real roots of / over an 
interval [x, y] is the difference of the variation of / at x and at y. 
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The Macaulay 2 commands numRealSturm and numPosRoots (and also numNegRoots) 
use this method to respectively compute the total number of real roots and the number 
of positive roots of a univariate polynomial. 

il8 : code ninnRealSturm 

ol8 = — code for numRealSturm: 
— realroots.m2: 161-165 
numRealSturm = f -> ( 

c : = SturmSequence f ; 
variations (signAtMinus Infinity \ c) 
- variations (signAt Infinity \ c) 



018 : Net 

119 : code numPosRoots 

019 = — code for numPosRoots: 

— realroots.m2: 170-174 
numPosRoots = f -> ( 

c := SturmSequence f ; 
variations (signAtZero \ c) 

- variations (signAt Infinity \ c) 

019 : Net 

These use the commands signAt* (f), which give the sign of f at *. (Here, * is one of 
Infinity, zero, or Minusinf inity. Also variations (c) computes the number of sign 
changes in the sequence c. 

120 : code variations 

020 = — code for variations: 

— realroots.m2: 187-195 
variations = c -> ( 

n := 0; 
last := 0; 

scan(c, X -> if X =!= then ( 

if last < and x > or last > 

and x < then n = n+1; 
last = x; 

)); 

n) 

020 : Net 

A more sophisticated method to compute the number of real roots which can also give 
information about their location uses the rank and signature of the symmetric trace form. 
Suppose I C k[X] is a zero-dimensional ideal and set A :— k[X]/I. For h e k[X], set 
Sh{f,g) '■= trace(m/,jg). It is an easy exercise that Sh is a symmetric bilinear form on A. 
The procedure traceForm(h) in realroots.m2 computes this trace form Sj^. 

121 : code traceForm 

021 = — code for traceForm: 

— realroots.m2: 200-208 
traceForm = h -> ( 

assert ( dim ring h == ) ; 

b := basis ring h; 

k := coef f IcientRing ring h; 

mm := substitute (contract (transpose b, h * b ** b) , k) ; 
tr := matrix {apply (first entries b, x -> 

trace regularRep x)}; 
adjoint (tr * mm, source tr, source tr) 
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021 : Net 

The value of this construction is the following theorem. 

Theorem 2.2.2 (|^, |18|). Suppose C M and I is a zero- dimensional ideal in k[xi, . . . , Xn] 
and consider V(/) C C". Then, for h G k[xi, ... , the signature a{Sh) and rank p{Sh) 
of the bilinear form Sh satisfy 

a{Sh) = #{a G V(/) n M" : /i(a) > 0} - #{a G V(/) n M" : /i(a) < 0} 
p{Sh) = #{a G V(/) : %) ^ 0} . 

That is, the rank of Sh counts roots in C" — V{h), and its signature counts the real 
roots weighted by the sign of h (which is —1, 0, or 1) at each root. The command 
traceFormSignature (h) in realroots.m2 returns the rank and signature of the trace 
form Sh- 

122 : code traceFormSignature 

022 = — code for traceFormSignature: 

— realroots.m2: 213-224 
traceFormSignature = h -> ( 

A := ring h; 
assert ( dim A == ) ; 
assert ( char A == ); 
S := QQ[Z] ; 

TrF := traceForm(h) S; 

IdZ := Z * id_ (S" (numgens source TrF)); 

f := detCTrF - IdZ) ; 

« "The trace form S_h with h = " « h « 

" has rank " << rank(TrF) << " and signature " << 
numPosRoots(f ) - numNegRoots (f ) << endl 

022 : Net 

The Macaulay 2 command numRealTrace(A) simply returns the number of real roots of 
/, given A = k[X]/I. 

123 : code numRealTrace 

023 = — code for numRealTrace: 

— realroots.m2: 229-237 
numRealTrace = A -> ( 

assert ( dim A == ); 
assert ( char A == ); 
S := QQ[Z] ; 

TrF := traceForm(l_A) ** S; 

IdZ := Z * id_ (S" (numgens source TrF)); 

f := detCTrF - IdZ) ; 

numPosRoots (f) -numNegRoots (f ) 

o23 : Net 

Example 2.2.3. We illustrate these methods on the following polynomial system. 

124 : R = QQ[x, y] ; 

125 : I = ideal (1 - x~2*y + 2*x*y~2, y - 2*x - x*y + x~2) ; 

025 : Ideal of R 

The ideal / has dimension zero and degree 5. 

126 : dim I, degree I 

026 = (0, 5) 



o26 : Sequence 
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We compare the two methods to compute the ehminant of x in the ring R/I. 

i27 : A = R/I; 

±28 : time g = eliminant (x, QQ [Z] ) 
— used 0.03 seconds 



4 3 2 
5Z +62 + Z - 2Z + 1 



o28 = Z 
o28 : QQ [Z] 



i29 : time g = charPolyCx, Z) 
— used 0.01 seconds 



o29 = Z 



4 3 2 
5Z + 6Z + Z 



2Z + 1 



029 : QQ [Z] 

The ehminant has 3 real roots, which we test in two different ways. 

i30 : numRealSturm(g) , numRealTrace (A) 

030 = (3, 3) 
o30 : Sequence 



We use Theorem p.2.2| to isolate these roots in the x, y-plane. 

131 : traceFormSignature(x*y) ; 

The trace form S_h with h = x*y has rank 5 and signature 3 

Thus all 3 real roots lie in the first and third quadrants (where xy > 0). We isolate these 
further. 

132 : traceFormSignature (x - 2); 

The trace form S_h with h = x - 2 has rank 5 and signature 1 
This shows that two roots lie in the first quadrant with x > 2 and one lies in the third. 
Finally, one of the roots lies in the triangle y > 0, x > 2, and x + y < 3. 



±33 : traceFormSignature (x + y - 3) ; 

The trace form S_h with h = x + y - 3 has rank 5 and signature -1 

X + y = 



Figure |I] shows these three roots (dots), as well as the lines 



3 and x = 2. 



X 



-1 



— 1 



X + y = 3 




Figure 1. Location of roots 



2.5. Homotopy methods. We describe symbolic-numeric homotopy continuation meth- 
ods for finding approximate complex solutions to a system of equations. These exploit 
the traditional principles of conservation of number and specialization from enumerative 
geometry. 

Suppose we seek the isolated solutions of a system F{X) = where F = (/i, . . . , fn) 
are polynomials in the variables X = {xi, . . . ,xn). First, a homotopy H{X,t) is found 
with the following properties: 
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1. H{X,1) = F{X). 

2. The isolated solutions of the start system H{X, 0) = are known. 

3. The system H{X,t) = defines finitely many (complex) curves, and each isolated 
solution of the original system F{X) = is connected to an isolated solution crj(O) 
of H{X, 0) = along one of these curves. 

Next, choose a generic smooth path 7(t) from to 1 in the complex plane. Lifting 7 to 
the curves H{X,t) = gives smooth paths ai(t) connecting each solution crj(O) of the 
start system to a solution of the original system. The path 7 must avoid the finitely many 
points in C over which the curves are singular or meet other components of the solution 
set H{X, t) = 0. 

Numerical path continuation is used to trace each path crj(t) from t = to t = 1. When 
there are fewer solutions to F{X) = than to H{X, 0) = 0, some paths will diverge or 
become singular as t —>■ 1, and it is expensive to trace such a path. The homotopy is 
optimal when this does not occur. 

When N = n and the are generic, set G{X) := {gi, . . . ,gn) with gi = (xj — l){xi — 
2) ■ ■ ■ {xi — di) where di := deg(/j). Then the Bezout homotopy 

H{X,t) := tF{X) + {l-t)G{X) 

is optimal. This homotopy furnishes an effective demonstration of the bound in Bezout's 
Theorem for the number of solutions to F{X) = 0. 

When the polynomial system is deficient, the Bezout homotopy is not optimal. When 
n > N (often the case in geometric examples), the Bezout homotopy does not apply. In 
either case, a different strategy is needed. Present optimal homotopies for such systems 
all exploit some structure of the systems they are designed to solve. The current state-of- 
the-art is described in |27[| . 



Example 2.2.4. The Grobner homotopy is an optimal homotopy that exploits a 
square-free initial ideal. Suppose our system has the form 

F := g^{X),... ,gUX), Ai(X), . . . , Arf(X) 

where gi{X), . . . , gm{X) form a Grobner basis for an ideal / with respect to a given term 
order -<, Ai, ... , A^ are linear forms with d = dim(V(/)), and we assume that the initial 
ideal in^J is square- free. This last, restrictive, hypothesis occurs for certain determinantal 
varieties. 

As in [^, Chapter 15], there exist polynomials gi{X,t) interpolating between gi{X) and 
their initial terms m^gi{X) 

giiX;l) = gi{X) and g,{X-{)) = m^gi{X) 

so that ((7i(X, t), . . . ,gm{X,t)) is a flat family with generic flbre isomorphic to / and 
special flbre in^J. The Grobner homotopy is 

H{X,t) := g^{X,t),... ,gm{X,t), Ai(X), . . . , A,(X). 

Since in_<J is square-free, V(in_<J) is a union of deg(J)-many coordinate (i-planes. We 
solve the start system by linear algebra. This conceptually simple homotopy is in general 
not efficient as it is typically overdetermined. 

3. Some enumerative geometry 



We use the tools we have developed to explore the enumerative geometric problems of 
cylinders meeting 5 general points and lines tangent to 4 spheres. 
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3.1. Cylinders meeting 5 points. A cylinder is the locus of points equidistant from 
a fixed line in M'^. The Grassmannian of lines in 3-space is 4-dimensional, which implies 
that the space of cylinders is 5-dimensional, and so we expect that 5 points in will 
determine finitely many cylinders. That is, there should be finitely many lines equidistant 
from 5 general points. The question is: How many cylinders/lines, and how many of them 
can be real? 



Bottema and Veldkamp P show there are 6 complex cylinders and Lichtblau [0 ob- 
serves that if the 5 points are the vertices of a bipyramid consisting of 2 regular tetrahedra 
sharing a common face, then all 6 will be real. We check this reality on a configuration 
with less symmetry (so the Shape Lemma holds). 

If the axial line has direction V and contains the point P (and hence has parameter- 
ization P + tY) , and if r is the squared radius, then the cylinder is the set of points X 
satisfying 



= r- 



2 



Expanding and clearing the denominator of HVp yields 

(5) = r||Vf + [V-(X-P)]2- ||X-Pf ||V||2. 

We consider cylinders containing the following 5 points, which form an asymmetric bipyra- 
mid. 

134 : Points = {{2, 2, }, {1, -2, 0}, {-3, 0, 0}, 

{0, 0, 5/2}, {0, 0, -3»; 

Suppose that P = (0, yn, 2/12) and V = (1,2/21,2/22)- 

135 : R = QQ[r, yll, yl2, y21, y22] ; 

136 : P = matrix{{0, yll, yl2»; 

1 3 

036 : Matrix R < R 

137 : V = matrix{{l, y21, y22»; 

1 3 

037 : Matrix R < R 

We construct the ideal given by evaluating the polynomial (^ at each of the five points. 

138 : Points = matrix Points ** R; 

5 3 

038 : Matrix R < R 

139 : I = ideal apply(0..4, i -> ( 

X := Points~{i}; 
r * (V * transpose V) + 
((X - P) * transpose V)"2) - 

((X - P) * transpose (X - P)) * (V * transpose V) 
) ! 

039 : Ideal of R 

This ideal has dimension and degree 6. 

140 : dim I , degree I 

040 = (0, 6) 
o40 : Sequence 

There are 6 real roots, and they correspond to real cylinders (with r > 0). 

141 : A = R/I; numPosRoots(charPoly(r, Z)) 
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o42 



3.2. Lines tangent to 4 spheres. We now ask for the lines liaving a fixed distance from 
4 general points. Equivalently, tliese are tlie lines mutually tangent to 4 spheres. Since 
the Grassmannian of lines is four-dimensional, we expect there to be only finitely many 



such lines. Macdonald, Pach, and Theobald [17] show that there are indeed 12 lines, and 



that all 12 may be real. This problem makes geometric sense over any field k not of 
characteristic 2, and the derivation of the number 12 is also valid for algebraically closed 
fields not of characteristic 2. 

A sphere in A;^ is V(g(l,x)), where g is a quadratic form on A;^ and x G k^. If our 
field does not have characteristic 2, then there is a symmetric 4x4 matrix M such that 
g(u) = uMu*. 

A line i having direction V and containing the point P is tangent to the sphere defined 
by q when the univariate polynomial in s 

g((l,P) + s(0,V)) = g(l,P) + 2s(l,P)M(0,V)* + s2g(0,V), 

has a double root. Thus its discriminant vanishes, giving the equation 

(6) ((1,P)M(0, V)*)^ - (1,P)M(1,P)* ■ (0, V)M(0, V)* = 0. 

The matrix M of the quadratic form q of the sphere with center (a, b, c) and squared 
radius r is constructed by Sphere (a,b,c,r). 

143 : Sphere = (a, b, c, r) -> ( 

matrix-[{a~2 + b~2 + c~2 - r ,-a ,-b ,-c }, 
{ -a ,1,0,0}, 
i -b ,0,1,0}, 
{ -c ,0,0,1}} 

); 

If a line i contains the point P = {0,yu,yi2) and i has direction V = (1, ?/2i, ^22), then 
tangentTo(M) is the equation for i to be tangent to the quadric uMu^ = determined 
by the matrix M. 

144 : R = QQ[yll, yl2, y21, y22] ; 

145 : tangentTo = (M) -> ( 

P := matrix{{l, 0, yll, yl2}}; 

V := matrix-C-CO , 1, y21, y22}}; 

(P * M * transpose V)~2 - 

(P * M * transpose p) * (V * M * transpose V) 

The ideal of lines having distance y/b from the four points (0, 0, 0), (4, 1, 1), (1, 4, 1), and 
(1, 1,4) has dimension zero and degree 12. 

146 : I = ideal (tangentTo (Sphere (0,0,0, 5) ) , 

tangentTo (Sphere (4,1,1,5)), 
tangentTo (Sphere (1,4, 1,5)) , 
tangentTo (Sphere (1,1, 4, 5))) ; 

046 : Ideal of R 

147 : dim I , degree I 

047 = (0, 12) 
o47 : Sequence 

Thus there are 12 lines whose distance from those 4 points is \/5. We check that all 12 
are real. 

148 : A = R/I; 
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149 : numRealSturm(eliminant(yll - yl2 + y21 + y22, QQCZ])) 
o49 = 12 

Since no eliminant given by a coordinate function satisfies the hypotheses of the Shape 
Lemma, we took the ehminant with respect to the hnear form yu — yi2 + ?/2i + 2/22- 

This example is an instance of Lemma 3 of |T^. These four points define a regular 
tetrahedron with volume V = 9 where each face has area A = and each edge has 

length e = -\/l8- That result guarantees that all 12 lines will be real when e/2 < r < 
A^/SV, which is the case above. 

4. Schubert calculus 

The classical Schubert calculus of enumerative geometry concerns linear subspaces hav- 
ing specified positions with respect to other, fixed subspaces. For instance, how many lines 
in meet four given lines? (See Example [4.4.2| .) More generally, let 1 < r < n and sup- 
pose that we are given general linear subspaces Li, . . . , of k"' with dim Lj = n—r+l—li. 
When h + ■ ■ ■ + Im = r{n — r), there will be a finite number d{r, n;li,... , Im) of r-planes 
in fc" which meet each Li non-trivially. This number may be computed using classical 
algorithms of Schubert and Fieri (see [ll5|). 

The condition on r-planes to meet a fixed (n—r +l—/)-plane non-trivially is called a 
(special) Schubert condition, and we call the data (r, n; /i, . . . , (special) Schubert data. 
The (special) Schubert calculus concerns this class of enumerative problems. We give two 
polynomial formulations of this special Schubert calculus, consider their solutions over M, 
and end with a question for fields of arbitrary characteristic. 

4.1. Equations for the Grassmannian. The ambient space for the Schubert calculus 
is the Grassmannian of r-planes in /c", denoted Gr,n- For H G Gr,n, the rth exterior 
product of the embedding H ^ k^ gives a line 

k ~ A'H — > ~ k(") . 

This induces the Fliicker embedding Gr,n ^ P^'-)"^. If H is the row space of an r by n 
matrix, also written H, then the Fliicker embedding sends H to its vector of (") maximal 
minors. Thus the r-subsets of {0, . . . , n—1}, Yr^n '■= subsets(n, r), index Fliicker coordi- 
nates of Gr,n- The Fliicker ideal of Gr,n is therefore the ideal of algebraic relations among 
the maximal minors of a generic r hj n matrix. 

We create the coordinate ring k[pa \ a E ^2^5] of P^ and the Fliicker ideal of G2,5. The 
Grassmannian Gr,n of r-dimensional subspaces of fc" is also the Grassmannian of r— 1- 
dimensional affine subspaces of P"~^. Macaulay 2 uses this alternative indexing scheme. 

150 : R = ZZ/lOl [apply (subsets (5, 2) , i -> p_i )] ; 

151 : I = Grassmannland , 4, R) 

o51 = Ideal (p p -p p +p p ,p •■■ 

{2, 3} {1, 4> {1, 3> {2, 4} {1, 2} {3, 4} {2, 3> • • • 

051 : Ideal of R 

This projective variety has dimension 6 and degree 5 

152 : dlm(Proj (R/D) , degree (I) 

052 = (6, 5) 
o52 : Sequence 
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This ideal has an important combinatorial structure |2^, Example 11.9]. We write each 
a G Yr,n as an increasing sequence a: ai < ■ ■ ■ < ar- Given a,/3 e Y^^^, consider the 
two-rowed array with a written above /3. We say a < (3 ii each column weakly increases. 
If we sort the columns of an array with rows a and /3, then the first row is the meet a A (3 
(greatest lower bound) and the second row the join a\J (3 (least upper bound) of a and 
(3. These definitions endow Y^.^ with the structure of a distributive lattice. Figure shows 

34 

\ 

24 

/ \ 
23 14 

\ X \ 
13 04 




02 



01 



Figure 2. Y2,5 

We give k\pa\ the degree reverse lexicographic order, where we first order the variables 
Pa by lexicographic order on their indices a. 

Theorem 4.4.1. The reduced Grohner basis of the Pliicker ideal with respect to this degree 
reverse lexicographic term order consists of quadratic polynomials 

g{a,(3) = Pa-Pp - Pav/3-PaA/3 + lowcr tcrms in -< , 

for each incomparable pair a,P in Yr^n, and all lower terms Xp^ ■ ps in g{a,P) satisfy 
7 < a A /5 and a\J [3 <5. 

The form of this Grobner basis implies that the standard monomials are the sortable 
monomials, those PaPp ■ ■ - p^ with a < (3 < ■ ■ ■ < 'j. Thus the Hilbert function of Gr,n 
may be expressed in terms of the combinatorics of Y,, „. For instance, the dimension of 
Gr,n is the rank of Y^^n, and its degree is the number of maximal chains. From Figure |^, 
these are 6 and 5 for Y2^5, confirming our previous calculations. 

Since the generators g{a,(3) are linearly independent, this Grobner basis is also a min- 
imal generating set for the ideal. The displayed generator in o51, 

P{2,3}P{1,4} - P{1,3}P{2,4} - P{1,2}P{3,4} , 

is ^(23,14), and corresponds to the underlined incomparable pair in Figure ^. Since 
there are 5 such incomparable pairs, the Grobner basis has 5 generators. As 62,5 has 
codimension 3, it is not a complete intersection. This shows how the general enumerative 
problem from the Schubert calculus gives rise to an overdetermined system of equations 
in this global formulation. 

The Grassmannian has a useful system of local coordinates given by Matr^n-r as follows 

(7) Y G Matr,n-r I ^ rOWSpaCe [Ir ■ Y] G Gr,n ■ 
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Let Lbea(n — r + 1 — /)-plane in k"^ which is the row space ofan — r + 1 — /byn 
matrix, also written L. Then L meets X G Gr^n non-trivially if 



maximal minors of 



L 
X 



0. 



Laplace expansion of each minor along the rows of X gives a linear equation in the Pliicker 
coordinates. In the local coordinates (substituting [J^ : Y] for X), we obtain multilinear 
equations of degree min{r, n — r}. These equations generate a prime ideal of codimension 
1. 

Suppose each /j = 1 in our enumerative problem. Then in the Pliicker coordinates, 
we have the Pliicker ideal of Gr,n together with r{n — r) linear equations, one for each 
(ra— r)-plane Lj. By Theorem |4.4.1| , the Pliicker ideal has a square-free initial ideal, and so 
the Grobner homotopy of Example |2.2.4| may be used to solve this enumerative problem. 



Example 4.4.2. 62,4 C has equation 

(8) P{1,2}P{0,3} - P{1,3}P{0,2} + P{2,3}P{0,1} = 0. 

The condition for H G 62,4 to meet a 2-plane L is the vanishing of 

(9) P{l,2}-^34 — P{l,3}-^24 + P{2,3}-^14 + P{l,4}-^23 — P{2,4}-^^13 + P{3,4}-^12 , 

where is the (z, j)th maximal minor of L. 

If li = ■ ■ ■ = I4 = 1, we have 5 equations in P^, one quadratic and 4 linear, and so by 
Bezout's Theorem there are two 2-planes in fc^ that meet 4 general 2-planes non-trivially. 
This means that there are 2 lines in P^ meeting 4 general lines. In local coordinates, 
becomes 

L34 - Luyn + Li3?/i2 - /^24y2i + -^^23^/22 + -^^i2(z/iil/22 - I/12I/21) • 
This polynomial has the form of the last specialization in Example |1.1.2| . 

4.2. Reality in the Schubert calculus. Like the other enumerative problems we have 
discussed, enumerative problems in the special Schubert calculus are fully real in that all 
solutions can be real [^. That is, given any Schubert data (r, n; /i, . . . , there exist 
subspaces Li, . . . , C M" such that each of the d{r, n-Ji, . . . , Z^) r-planes that meet 
each Li are themselves real. 

This result gives some idea of which choices of the Li give all r-planes real. Let 
7 be a fixed rational normal curve in M". Then the Lj are linear subspaces osculat- 
ing 7. More concretely, suppose that 7 is the standard rational normal curve, 7(5) = 
(1, s, s^, . . . ,s"'~^). Then the z-plane Li(s) := (7(5), 7'(s), . . . ,7'^*~^^(s)) osculating 7 at 
7(5) is the row space of the matrix given by oscPlane(i, n, s). 

153 : oscPlane = (i, n, s) -> ( 

gamma := matrix {toList apply(l..ii, i -> s"(i-l))}; 
L := gamma; 
j := 0; 

while j < i-1 do (gsunma = diff(s, gamma); 
L = L I I gamma; 

j = 

L); 

154 : QQ[s]; oscPlane (3, 6, s) 

o55 = I 1 s s2 s3 s4 s5 I 
I 1 2s 3s2 4s3 5s4 I 
10 2 6s 12s2 20s3 I 



3 6 
o55 : Matrix QQ [s] < QQ [s] 



16 



FRANK SOTTILE 



(In o55, the exponents of s are displayed in line: is written s2. Macaulay 2 uses this 
notational convention to display matrices efficiently.) 

Theorem 4.4.3 (|^). For any Schubert data (r, n;li, . . . , there exist real numbers 
Si,S2, ■ ■ ■ ,Sm such that there are d{r, n;li,... , 1^) r -planes that meet each osculating plane 
Li{si), and all are real. 

The inspiration for looking at subspaces osculating the rational normal curve to study 
real enumerative geometry for the Schubert calculus is the following very interesting con- 
jecture of Boris Shapiro and Michael Shapiro, or more accurately, extensive computer 



experimentation based upon their conjecture [T^, ^ . 

Shapiros's Conjecture. For any Schubert data (r, n]li, . . . , 1^) and for all real num- 
bers si,S2,... ,Sm there are d{r,n]li,... Jm) r -planes that meet each osculating plane 
Li{si), and all are real. 

In addition to Theorem ^.4.31 , (which replaces the quantifier for all by there exist), the 
strongest evidence for this Conjecture is the following result of Eremenko and Gabrielov p!0| 



Theorem 4.4.4. Shapiros's Conjecture is true when either r or n — r is 2. 

We test an example of this conjecture for the Schubert data (3, 6; 1^, 2^), (where is a 
repeated b times). The algorithms of the Schubert calculus predict that d{3, 6; 1^, 2^^) = 6. 
The function spSchub(r, L, P) computes the ideal of r-planes meeting the row space of 
L in the Pliicker coordinates Pa- 

156 : spSchub = (r, L, P) -> ( 

I := ideal apply (subsets (numgens source L, 

r + numgens target L) , S -> 
fold ((sum, U) -> sum + 
fold((term,i) -> term*(-l)~i, P_(S_U) * det ( 
submatrix(L, sort toList (set (S) - set(S_U)))), U) , 
0, subsets (#S, r)))) ; 

We are working in the Grassmannian of 3-planes in C^. 

157 : R = QQ[apply(subsets(6,3) , i -> p_i )]; 

The ideal / consists of the special Schubert conditions for the 3-planes to meet the 3- 
planes osculating the rational normal curve at the points 1, 2, and 3, and to also meet the 
2-planes osculating at 4, 5, and 6, together with the Pliicker ideal Grassmannian (2 , 5, 
R). Since this is a 1-dimensional homogeneous ideal, we add the linear form p_{0,l,5} 
- 1 to make the ideal zero-dimensional. As before, Grassmannian (2, 5, R) creates the 
Pliicker ideal of G^^. 

158 : 1 = fold((J, 1) -> J + 

spSchub(3, substitute (oscPlane (3 , 6, s) , {s=> 1+1}), p) + 

spSchub(3, substitute (oscPlane (2 , 6, s) , {s=> 4+1}), p) , 
Grassmannian (2, 5, R) , {0,1,2}) + 
Ideal (p_{0,l,5} - 1) ; 

058 : Ideal of R 

This has dimension and degree 6, in agreement with the Schubert calculus. 

159 : dim 1, degree 1 

059 = (0, 6) 
o59 : Sequence 

As expected, all roots are real. 

160 : A = R/1; numRealSturm(ellmlnant(p_{2,3,4}, QQ[Z])) 



o61 = 6 
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There have been many checked instances of this conjecture p3| , and it has some 



geometrically interesting generalizations 24 



The question remains for which numbers < d < d{r, n;li, . . . , Im) do there exist real 
planes Li with d{r, n;li, . . . , /„) r-planes meeting each Li, and exactly c? of them are real. 
Besides Theorem [4.4.3| and the obvious parity condition, nothing is known in general. 



In every known case, every possibility occurs — which is not the case in all enumerative 
problems, even those that are fully realQ. Settling this (for d = 0) has implications for 
linear systems theory ||TU|. 

4.3. Transversality in the Schubert calculus. A basic principle of the classical Schu- 
bert calculus is that the intersection number d{r, n;li,... , Im) has enumerative significance- 
that is, for general linear subspaces Lj, all solutions appear with multiplicity 1. This 
basic principle is not known to hold in general. For fields of characteristic zero, Kleiman's 
Transversality Theorem establishes this principle. When r or n—r is 2, then Theo- 



rem E of ||2^ establishes this principle in arbitrary characteristic. We conjecture that this 
principle holds in general; that is, for arbitrary infinite fields and any Schubert data, if 
the planes Li are in general position, then the resulting zero- dimensional ideal is radical. 

We test this conjecture on the enumerative problem of Section |4.2| , which is not covered 
by Theorem E of [^]. The function testTransverse (F) tests transversality for this 



enumerative problem, for a given field F. It does this by first computing the ideal of the 
enumerative problem using random planes Lj. 

162 : randL = (R, n, r, 1) -> 

matrix table (n-r+l-l , n, (i, j) -> randomCO, R) ) ; 

and the Pliicker ideal of the Grassmannian Gs^ Grassmannian(2 , 5, R).) Then it adds 
a random (inhomogeneous) linear relation 1 + randomd, R) to make the ideal zero- 
dimensional for generic Lj. When this ideal is zero dimensional and has degree 6 (the 
expected degree) , it computes the characteristic polynomial g of a generic linear form. If 
g has no multiple roots, 1 == gcd(g, diff (Z, g)), then the Shape Lemma guarantees 
that the ideal was radical. testTransverse exits either when it computes a radical ideal, 
or after limit iterations (which is set to 5 for these examples), and prints the return 
status. 

163 : testTransverse = F -> ( 

R := F [apply (subsets (6, 3), i -> q_i )] ; 
continue := true; 
j := 0; 
limit := 5; 

while continue and (j < limit) do ( 

j = j + 1; 

I := fold((J, i) -> J + 

spSchub(3, randLCR, 6, 3, 1), q) + 

spSchub(3, randLCR, 6, 3, 2), q) , 

Grassmannian (2 , 5, R) + 

ideal (1 + randomCl, R) ) , 

{0, 1, 2»; 
if (dim I == 0) and (degree I == 6) then ( 
lin := promote (randomd , R) , (R/I)); 
g := charPoly (lin, Z) ; 

continue = not(l == gcd(g, diff(Z, g) ) ) ; 
) ) > 

if continue then << "Failed for the prime " << char F << 

" with " « j « " iterations" « endl; 
if not continue then << "Succeeded for the prime " << 
char F « " in " « j « " iteration(s) " « endl; 



^For example, of the 12 rational plane cubics containing 8 real points in P^, either 8, 10 or 12 can be 
real, and there are 8 points with all 12 real |@, Proposition 4.7.3]. 
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); 

Since 5 iterations do not show transversality for F2, 

164 : testTransverse(ZZ/2) ; 

Failed for the prime 2 with 5 iterations 

we can test transversality in characteristic 2 using the field with four elements, F4 = GF 
4. 

165 : testTransverse (GF 4); 

Succeeded for the prime 2 in 3 iteration(s) 

We do find transversality for F7. 

166 : testTransverse (ZZ/7) ; 

Succeeded for the prime 7 in 2 iteration(s) 

We have tested transversality for all primes less than 100 in every enumerative problem 
involving Schubert conditions on 3-planes in k^. These include the problem above as well 
as the problem of 42 3-planes meeting 9 general 3-planes.0 

5. The 12 lines: reprise 

The enumerative problems of Section ^ were formulated in local coordinates (0) for the 
Grassmannian of lines in (Grassmannian of 2-dimensional subspaces in fc^). When we 

we 



formulate the problem of Section in the global Pliicker coordinates of Section |4.1| , 
find some interesting phenomena. We also consider some related enumerative problems. 

5.1. Global formulation. A quadratic form g on a vector space V over a field k not 
of characteristic 2 is given by g(u) = ((y9(u),u), where (p: V ^ V* is a symmetric linear 
map, that is (y9(u), v) = {{p{'v),u). Here, V* is the linear dual of V and (■ , ■ ) is the 
pairing V ^V* k. The map induces a quadratic form A^q on the rth exterior power 
A^'V of V through the symmetric map A^'p: A"^ V A^'V* = (A^'V)*. The action of A'^V* 
on A^V is given by 

(10) (xi A X2 A • ■ ■ A x^, A ys A ■ ■ ■ A y^) = det |(xi, y-)| , 

where Xj G V^* and jj G V. 

When we fix isomorphisms ~ A;" ~ V^*, the map p is given by a symmetric n x n 
matrix M as in Section |3.2| . Suppose r = 2. Then for u, v G /c", 

2 , , , , r uMu* uMv* 
A g(u A v) = det 

which is Equation (B) of Section ^72. 



vMu* vMv* 



Proposition 5.5.1. A line i is tangent to a quadric V{q) in ^ if and only if its Pliicker 
coordinate A^l G p(2)~^ lies on the quadric V(A^g). 

Thus the Pliicker coordinates for the set of lines tangent to 4 general quadrics in P^ 
satisfy 5 quadratic equations: The single Pliicker relation together with one quadratic 
equation for each quadric. Thus we expect the Bezout number of 2^ = 32 such lines. We 
check this. 

The procedure randomSymmetricMatrix(R, n) generates a random symmetric n x n 
matrix with entries in the base ring of R. 

i67 : randomSymmetricMatrix = (R, n) -> ( 
entries := new MutableHashTable ; 
scan(0..n-l, i -> scan(i..n-l, j -> 



2 



After this was written, we discovered an elementary proof of transversality for the enumerative prob- 



lems d{r,n] f^" '"'), where the conditions are all codimension 1 |25 
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entries* (i, j) = random(0, R) ) ) ; 
matrix table (n, n, (i, j) -> if i > j then 

entries#(j, i) else entries#(i, j)) 

) > 

The procedure tangentEquationCr , R, M) gives the equation in Pliicker coordinates for 

a point in p(')^^ to be isotropic with respect to the bihnear form A''M (R is assumed to 

be the coordinate ring of p(')~^). This is the equation for an r-plane to be tangent to the 
quadric associated to M. 

168 : tangentEquation = (r, R, M) -> ( 

g := matrix {gens(R)}; 

(entries (g * exteriorPower (r , M) * transpose g))_0_0 

We construct the ideal of lines tangent to 4 general quadrics in P'^. 

169 : R = QQ[apply(subsets(4, 2), 1 -> p_i )] ; 

170 : I = GrassmannianCl , 3, R) + ideal apply(0..3, 1 -> 

tangentEquation(2 , R, rcindomSymmetrlcMatrlxCR, 4))); 

070 : Ideal of R 

As expected, this ideal has dimension and degree 32. 

171 : dim Proj(R/l), degree I 

071 = (0, 32) 
o71 : Sequence 



5.2. Lines tangent to 4 spheres. That calculation raises the following question: In 



Section 3.2, why did we obtain only 12 lines tangent to 4 spheres? To investigate this, we 



generate the global ideal of lines tangent to the spheres of Section ^ 

172 : 1 = GrassmannianCl, 3, R) + 

ideal (tangentEquation(2 , R, Sphere (0 , , , 5) ) , 
tangentEquation(2 , R, Sphere (4, 1 , 1 , 5) ) , 
tangentEquation(2 , R, Sphere (1 ,4, 1 , 5) ) , 
tangentEquation(2 , R, Sphere (1 , 1 ,4, 5) ) ) ; 

072 : Ideal of R 

We compute the dimension and degree of V(/). 

173 : dim Proj (R/1) , degree 1 

073 = (1, 4) 

073 : Sequence 

The ideal is not zero dimensional; there is an extraneous one-dimensional component of 
zeroes with degree 4. Since we found 12 lines in Section ^]2| using the local coordinates (^, 
the extraneous component must lie in the complement of that coordinate patch, which is 
defined by the vanishing of the first Pliicker coordinate, P{o,i}- We saturate / by P{o,i} to 
obtain the desired lines. 

174 : Lines = saturated, ideal (p_{0,l})); 

074 : Ideal of R 

This ideal does have dimension and degree 12, so we have recovered the zeroes of 
Section 



175 : dim Proj (R/Llnes) , degree (Lines) 
o75 = (0, 12) 



o75 : Sequence 
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We investigate the rest of the zeroes, which we obtain by taking the ideal quotient of / 
and the ideal of lines. As computed above, this has dimension 1 and degree 4. 

176 : Junk = I : Lines; 

076 : Ideal of R 

177 : dim Pro j (R/ Junk) , degree Junk 

077 = (1, 4) 

077 : Sequence 

We find the support of this extraneous component by taking its radical. 

178 : radical (Junk) 

2 2 2 

078 = Ideal (p ,p ,p ,p +p +p ) 

{0, 3} {0, 2> {0, 1> {1, 2> {1, 3> {2, 3> 

o78 : Ideal of R 

From this, we see that the extraneous component is supported on an imaginary conic in 
the of lines at infinity. 

To understand the geometry behind this computation, observe that the sphere with 
radius r and center (a, 6, c) has homogeneous equation 

(x — waY + iy — whY + {z — wcf = r'^w'^ . 

At infinity, w = 0, this has equation 

x'^ + y'^ + z'^ = . 

The extraneous component is supported on the set of tangent lines to this imaginary conic. 
Aluffi and Fulton studied this problem, using geometry to identify the extraneous ideal 
and the excess intersection formula to obtain the answer of 12. Their techniques show 
that there will be 12 isolated lines tangent to 4 quadrics which have a smooth conic in 
common. 

When the quadrics are spheres, the conic is the imaginary conic at infinity. Fulton 
asked the following question: Can all 12 lines be real if the (real) four quadrics share a 
real conic? We answer his question in the affirmative in the next section. 

5.3. Lines tangent to real quadrics sharing a real conic. We consider four quadrics 
in Pj^ sharing a non-singular conic, which we will take to be at infinity so that we may use 
local coordinates for 62,4 in our computations. The variety V(g) C P| of a nondegenerate 
quadratic form q is determined up to isomorphism by the absolute value of the signature 
a of the associated bilinear form. Thus there are three possibilities, 0, 2, or 4, for \a\. 

When \a\ = 4, the real quadric V(g) is empty. The associated symmetric matrix M is 
conjugate to the identity matrix, so A^M is also conjugate to the identity matrix. Hence 
V(A^g) contains no real points. Thus we need not consider quadrics with \a\ = 4. 

When \a\ = 2, we have V{q) ^ S"^, the 2-sphere. If the conic at infinity is imaginary, 
then V{q) C M.^ is an ellipsoid. If the conic at infinity is real, then V{q) C M'^ is a 
hyperboloid of two sheets. When a = 0, we have V(g) ~ 5*^ x 5*^ a torus. In this case, 
V{q) C M'^ is a hyperboloid of one sheet and the conic at infinity is real. 

Thus either we have 4 ellipsoids sharing an imaginary conic at infinity, which we studied 
in Section |3.2| ; or else we have four hyperboloids sharing a real conic at infinity, and there 
are five possible combinations of hyperboloids of one or two sheets sharing a real conic at 
infinity. This gives six topologically distinct possibilities in all. 
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Theorem 5.5.2. For each of the six topologically distinct possibilities of four real quadrics 
sharing a smooth conic at infinity, there exist four quadrics having the property that each 
of the 12 lines in simultaneously tangent to the four quadrics is real. 



Proof. By the computation in Section p.2| , we need only check the five possibihties for 



hyperboloids. We fix the conic at infinity to be x + y — z 
hyperboloid of two sheets containing this conic has equation in 



0. Then the general 



'111 



(x - af + {y- hf -{z-cf + r 



0. 



(with r > 0). The command Two(a,b,c,r) generates the associated symmetric matrix. 

179 : Two = (a, b, c, r) -> ( 

matrlx{{a~2 + b~2 - c~2 + r ,-a ,-b , c }, 
{ -a ,1,0,0}, 

{ -b ,0,1,0}, 

{ c ,0,0,-1}} 

); 

The general hyperboloid of one sheet containing the conic a;^ + — 2;^ = at infinity has 
equation in 

(12) {x-af + {y-hf -{z-cf -r = 

(with r > 0). The command One(a,b,c,r) generates the associated symmetric matrix. 

180 : One = (a, b, c, r) -> ( 

matrlx{{a~2 + b~2 - c"2 
{ -a 
{ -b 
{ c 

); 

We consider i quadrics of two sheets ( pT]) and A — i quadrics of one sheet (|1^). For each 
of these cases, the table below displays four 4-tuples of data (a, 6, c, r) which give 12 real 
lines. (The data for the hyperboloids of one sheet are listed first.) 



r ,-a ,-b , c }, 
,1,0,0}, 
,0,1 
,0,0 



, }, 
,-1 }} 



i 


Data 





(5,3,3,16), 


(5,- 


-4,2,1), (- 


-3,-1,1,1), 


(2,-7,0,1) 


1 


(3,-2,-3,6), 


(-3,- 


-7,-6,7), ( 


-6,3, -5,2), 


(1,6,-2,5) 


2 


(6,4,6,4), 


(-1 


,3,3,6), (- 


-7, -2,3,3), 


(-6,7,-2,5) 


3 


(-1,-4,-1,1), 


(-3, 


3,-1,1), 


(-7,6,2,9), 


(5,6,-1,12) 


4 


(5,2,-1,25), 


(6,- 


-6,2,25), 


:-7, 1,6,1), 


(3,1,0,1) 



We test each of these, using the formulation in local coordinates of Section ^ 
181 : R = QQ[yll, yl2, y21, y22] ; 
I 



182 



Ideal (tangentTo(Dne( 5, 3, 3,16)), 
tangentTo(One( 5,-4, 2, 1)), 
tangentTo (One (-3,-1, 1, 1)), 
tangentTo(One( 2,-7, 0, 1))) 



082 : Ideal of R 

183 : numRealSturm(charPoly (promote (y22, R/I) , Z)) 

083 = 12 

184 : I = Ideal (tangentTo (One ( 3,-2,-3, 6)), 

tangentTo (One (-3, -7, -6, 7)), 
tangentTo (One (-6 , 3,-5, 2)), 
tangentTo (Two ( 1, 6,-2, 5))); 



o84 



Ideal of R 
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185 : numRealSturm(charPoly (promote (y22, R/I) , Z)) 

085 = 12 

186 : I = ideal (tangentTo (One ( 6, 4, 6, 4)), 

tangentTo (One (-1 , 3, 3, 6)), 
tangentTo (Two (-7, -2, 3, 3)), 
tangentTo (Two (-6 , 7,-2, 5))); 

086 : Ideal of R 

187 : numRealSturm(charPoly (promote (y22, R/1) , Z)) 

087 = 12 

188 : I = Ideal (tangentTo(0ne(-l,-4,-l, 1)), 

tangentTo (Two (-3, 3,-1, 1)), 
tangentTo (Two (-7, 6, 2, 9)), 
tangentTo (Two ( 5, 6,-1,12))); 

088 : Ideal of R 

189 : numRealSturm(charPoly (promote (y22, R/I), Z)) 

089 = 12 

190 : I = Ideal (tangentTo (Two ( 5, 2,-1,25)), 

tangentTo (Two ( 6,-6, 2,25)), 
tangentTo (Two (-7, 1, 6, 1)), 
tangentTo (Two ( 3, 1, 0, 1))); 

090 : Ideal of R 

191 : numRealSturm(charPoly (promote (y22, R/I), Z)) 

091 = 12 

□ 

In each of these enumerative problems, we have checked that every possible number of 
real solutions (0, 2, 4, 6, 8, 10, or 12) can occur. 

5.4. Generalization to higher dimensions. We consider lines tangent to quadrics in 
higher dimensions. First, we reinterpret the action of A'^V^* on A'^V described in (|TU|) 
as follows. The vectors xi, . . . , and yi, . . . , yr define maps a: —>■ V* and (3: ^ 
V. The matrix [(x,, yj)] is the matrix of the bilinear form on k^ given by (u, v) : = 
(q;(u), /9(v)). Thus (0) vanishes when the bilinear form ( ■ , ■ ) on fc^ is degenerate. 

Now suppose that we have a quadratic form qonV given by a symmetric map ip: V ^ 
V*. This induces a quadratic form and hence a quadric on any r-plane H in V (with 
H (t V(g)). This induced quadric is singular when H is tangent to V(g). Since a quadratic 
form is degenerate only when the associated projective quadric is singular, we see that H 
is tangent to the quadric V(g) if and only if {A^{p{A^H), A^H) = 0. (This includes the 
case H C V{q)-) We summarize this argument. 

Theorem 5.5.3. Let (f: V V* be a linear map with resulting bilinear form {f{u), v). 
Then the locus of r -planes in V for which the restriction of this form is degenerate is the 
set of r -planes H whose Pliicker coordinates are isotropic, {A^(f{A'^H), A^H) = 0, with 
respect to the induced form on A'^V. 

When (f is symmetric, this is the locus of r -planes tangent to the associated quadric in 
f{V). 

We explore the problem of lines tangent to quadrics in P". From the calculations of 
Section |5]5, we do not expect this to be interesting if the quadrics are general. (This is 
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borne out for P"^: we find 320 lines in tangent to 6 general quadrics. This is the Bezout 
number, as deg 62,5 = 5 and the condition to be tangent to a quadric has degree 2.) This 
problem is interesting if the quadrics in P" share a quadric in a P"~^. We propose studying 
such enumerative problems, both determining the number of solutions for general such 
quadrics, and investigating whether or not it is possible to have all solutions be real. 

We use Macaulay 2 to compute the expected number of solutions to this problem when 
r — 2 and n — A. We first define some functions for this computation, which will involve 
counting the degree of the ideal of lines in P^ tangent to 6 general spheres. Here, X 
gives local coordinates for the Grassmannian, M is a symmetric matrix, tanQuad gives 
the equation in X for the lines tangent to the quadric given by M. 

192 : tanQuad = (M, X) -> ( 

u := X-{0}; 
V := X-{1}; 

(u * M * transpose v)~2 - 

(u * M * transpose u) * (v * M * transpose v) 
) ! 

nSphere gives the matrix M for a sphere with center V and squared radius r, and V and 
r give random data for a sphere. 

193 : nSphere = (V, r) -> 

(matrix {{r + V * transpose V}} I I transpose V ) I 
(VII id_((ring r)-n) 
) > 

194 : V = () -> matrix tabled, n, (i,j) -> random(0, R)); 

195 : r = -> random(0, R) ; 

We construct the ambient ring, local coordinates, and the ideal of the enumerative problem 
of lines in P^ tangent to 6 random spheres. 

196 : n = 4; 

197 : R = ZZ/1009[flatten(table(2, n-1, (i,j) -> z_(i,j)))]; 

198 : X = 1 I matrix table (2, n-1, (i,j) -> z_(i,j)) 

o98 = I 1 z_(0,0) z_(0,l) z_(0,2) j 
I 1 z_(l,0) z_(l,l) z_(l,2) I 

2 5 

098 : Matrix R < R 

199 : I = ideal (apply (1 .. (2*n-2) , 

i -> tanC!uad(nSphere(V() , r()), X))); 

099 : Ideal of R 

We find there are 24 lines in P^ tangent to 6 general spheres. 

ilOO : dim I, degree I 

0100 = (0, 24) 
olOO : Sequence 

The expected numbers of solutions we have obtained in this way are displayed in the table 
below. The numbers in boldface are those which are proven. 



n 


2 


3 


4 


5 


6 


# expected 


4 


12 


24 


48 


96 
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